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We consider the scattering of two color dipoles (e.g., heavy quarkonium states) at low 
energy - a QCD analog of Van der Waals interaction. Even though the couplings of the 
ON ' dipoles to the gluon field can be described in perturbation theory, which leads to the potential 

ON ' proportional to (iV<? — 1)/R 7 , at large distances R the interaction becomes totally non- 

on : perturbative. Low-energy QCD theorems are used to evaluate the leading long-distance 

contribution ~ (iV| - 1)/(1UV C - 2N f ) 2 R~ 5/2 exp(~2^R) (fi is the Goldstone boson mass), 
(3JT)' which is shown to arise from the correlated two-boson exchange. The sum rule which relates 

the overall strength of the interaction to the energy density of QCD vacuum is derived. 
Surprisingly, we find that when the size of the dipoles shrinks to zero (the heavy quark 
limit in the case of quarkonia), the non-perturbative part of the interaction vanishes more 
slowly than the perturbative part as a consequence of scale anomaly. As an application, we 
evaluate elastic irJ/if) and nJ/ip inp' cross sections. 
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The interaction between small color dipolesQ provides an interesting theoretical laboratory for the studies 
of QCD and its applications in nuclear physics. Indeed, the asymptotic freedom dictates that the coupling of 
strong interactions becomes weak at short distances, and since the small size of dipoles introduces a natural 
infrared cut-off, one can hope that their interactions can be systematically treated in perturbation theory 
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D , One could therefore expect that at low energy the interaction between the dipoles in SU(N) gauge theory 

• would be of Van der Waals type: 



v peit (R) ~ ~9\N 2 - 1) (1) 

where n = 6 in the original Van der Waals potential, and g is the gluon coupling evaluated at the scale 
of quarkonium size. Indeed, this behavior was established by Appelquist and Fischler jj]], who studied the 
interactions of static color dipoles described by Wilson loops. These authors also explored the breakdown of 
the perturbative expansion in the static potential [||, and pointed to the possibility that retardation effects 
can modify the 1/i? 6 dependence once the spatial motion of the quarks is considered. In this paper, we 
take this effect into account and argue that, in the limit of the small size of the dipoles, the potential (|l|) 
is actually of Casimir-Polder type, with n = 7. On the other hand, gluons cannot propagate at large 
distances, where the dominant degrees of freedom are the lightest hadronic states. In the chiral limit, the 
theory with spontaneously broken SUi,(Nf) x SUa(Nf) symmetry contains (iVj — 1) Goldstone bosons, and 
the number of flavors Nf should effectively replace the number of colors in the coefficient of Eq. (|l|) at large 
distances: 

V cMral (R) ~ -{N} - 1) (2) 



*Present address: Institute of Physics, University of Tokyo, Komaba, Meguro, Tokyo 153-8902. 
1 Small color dipoles can be realized in the real word as heavy quarkonium states or as virtual quark-antiquark pairs 
in deep inelastic scattering. 
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In the real world where the masses [i of Goldstone bosons are not equal to zero, instead of Eq. (|^) at large 
distances one expects to find the potential of Yukawa form 
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V(R)~-(Nj-l) — — . (3) 

(We will show that the actual form of the long-distance potential is different from (||) - see Eq. ([n]).) How 
does the transition between the behavior at short and long distances occur? Can one explicitly, from the 
first principles, evaluate the long-distance potential? 

In this paper we address these questions, and argue that the interaction between small color dipoles 
(heavy quarkonium states in our example) at large distances can be reliably evaluated. Our analysis is 
based on the following two properties of QCD: 1) the scale invariance which is present at the tree level in 
QCD with massless quarks is broken by interactions; this is reflected in the non-zero divergence of scale 
current, and hence non- vanishing trace of the energy- momentum tensor [§U|]; 2) the chiral symmetry is 
broken spontaneously, which implies the existence of Goldstone bosons; being the lightest of all hadrons, 
they are the relevant degrees of freedom at large distances. 

These two properties of QCD are beautifully linked by the low-energy theorem derived by Voloshin and 
Zakharov flO|| , which we discuss below. The first of these properties was previously exploited to derive 
the low-energy amplitude of quarkonium-nucleon scattering jn],[l2| (for recent work, see |l3],[F4|] ) . Van der 
Waals interactions of quarkonium with nucleons and nuclei were discussed in Refs. |15|-1S|. For applications 
to the low-energy quarkonium-pion scattering and the structure of quarkonium, see Jl9-22|. Quarkonium 
dissociation cross sections in interactions with light hadrons were evaluated in Refs. [p| p3p4[| . Some of the 
results of this study were previously reported in Ref. j25[ . 

The picture which emerges from our approach is the following. The heavy "onia" couple perturbatively to 
the gluon field; at small distances, the entire interaction can be evaluated perturbatively. At larger distances, 
however, the interaction becomes grossly modified by the coupling to pion fields, which is fixed by low-energy 
theorems. The dominance of the non-perturbative interaction at large distances in this case will be shown to 
be a consequence of the finite energy density of QCD vacuum f| We also find that when the size of the dipoles 
shrinks to zero (which is what happens in the heavy quark limit with quarkonia), the non-perturbative part 
of the interaction vanishes more slowly than the perturbative part - in other words, the interaction between 
very small dipoles becomes totally non-perturbative! This surprising result will be shown to be a natural 
consequence of scale anomaly in QCD. 

In this paper, we limit ourselves to the interaction at small energies; however we hope that some of 
our results may be extended to the case of dipole scattering at high energies |P6|, where the broken chiral 
symmetry can also play a substantial role, as discussed by Anselm and Gribov ||27|. 

The paper is organized as follows. In Section [n], we give the general expression for the scattering amplitude 
of two color dipoles in the framework of the Operator Product Expansion (OPE), introduce the spectral 
representation method for the evaluation of this amplitude, and use this method to r e-derive the perturbative 
expression || for the low-energy scattering amplitude (or potential). In Section |n| we discuss the scattering 
amplitude of color dipoles beyond the perturbation theory, derive the leading long-distance behavior of the 
potential, and discuss the relative strength of perturbative and non-perturbative contributions. In Section [TV , 
we evaluate the potential acting between two J/t/Vs. In Section [v| we use the low energy theorems [ p~0|j20| pi ] 
to derive the sum rule relating the strength of the potential to the energy density of QCD vacuum. In 



Section VI, we evaluate the cross sections of J ftp interactio ns w ith pions, relevant for the problem of J/tp 



suppression in heavy ion collisions p2j,p3| . The final Section VII is devoted to summary and discussion 



II. INTERACTION OF COLOR DIPOLES IN PERTURBATION THEORY 

The small size of the heavy quarkonium $ allows us to expand the amplitude of its interaction with 
hadrons (h, h') at low energy in the form of multipole, or operator product, expansion 

M = Y^*W\OM\h), (4) 



2 This picture was foreseen by Bjorkcri 
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where Oi(x) are the gauge-invariant local operators and Ci are the Wilson coefficients (polarizabilities) which 
reflect the structure of the quarkonium; the energy of the hadrons is assumed to be small compared to the 
binding energy of the quarkonium, eo. The factorization scale in this formula can be chosen at eo. At small 
energies, the leading operator in Eq. (Q) is the square of the chromo-electric field (l/2)g 2 E a2 (0) this 
is the leading twist-two operator expressible in terms of gluon fields. Other twist-two operators contain 
covariant derivatives leading to the powers of the ratio of the energy transfer to the binding energy and are 
therefore suppressed at small energy; the series of these local operators can be summed up into a double 
dipole form H: 
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tr[^(0)H£T^(0)] = ^(tr r • E(0) — — — -r • E(0) ) , (5) 



H a + e + iD° 



where the Wilson coefficients are explicitly given by the expectation values over the singlet state cj>(r) with 
the binding energy e; H a (r) is the effective Hamiltonian describing the intermediate, SU(N) color-adjoint 
quark-antiquark state; D° is the covariant derivative acting on E and the trace over the color indices of 
gluon operators ensures that Eq. (j^) is gauge-invariant. In the heavy quark limit, <j)(r) can be approximated 
by the Coulomb wave function. 

Using this multipole representation, one can write down the amplitude of the scattering of two small color 
dipoles at low energies (in the Born approximation) in the following form gj: 

/oo 
dt (OlT^CiO^XX^COfllO). (6) 
j 

Keeping only the leading operators, we can rewrite Eq. (^|) in a simple form, 

d 2 Jl) / dt (0|T -g 2 E a • E a (i, R) -g 2 E b • E b (0)|0), (7) 
where d 2 is the corresponding Wilson coefficient defined by 

*- = 4<^?V ri i^ (8) 

e 3iV H a + e 

from which we have explicitly factored out the dependence on the quarkonium Bohr radius ao and the 
Rydberg energy eo- The factors ao and 1/eo represent the characteristic dimension and fluctuation time of 
the color dipole, respectively. The approximation used in deriving Eq. (^) is justified when the gluon fields 
change slowly compared to 1/eo- The Wilson coefficients (^) were computed for S || and P [|l9| states of 
quarkonium in the large TV limit. 

In physical terms, the structure of Eq. is transparent: it describes elastic scattering of two dipoles 
which act on each other by chromo-electric dipole fields; color neutrality permits only the square of dipole 
interaction. 

The amplitude (Q) was evaluated before [5j in perturbative QCD using functional methods. For our 
purposes, however, it is convenient to use a spectral representation approach^. As a first application of this 
approach, we will reproduce the result of || by a different, and perhaps more simple, method. 

First, it is convenient to express <? 2 E a2 in terms of the gluon field strength tensor fl29|| : 

ff 2 E a2 _ ^( E a2 _ B a2) + 9_ ^ a 2 + R a2) 

= --/GlpG™? + g 2 (-G a Qa G a a + \g m G a ap G^) = *?-9* + g 2 9^ , (9) 



where 



nil _ m 



3 The use of dispersion theory in electrodynamics for the interaction between neutral atoms was pioneered by Feinberg 
and Sucher 1341 . 
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Note that 6£ is the trace of the energy-momentum tensor of QCD in the chiral limit, and as a consequence 
of decoupling theorem |35| the f3 function in Eq. ( |Io|) docs not contain the contribution of heavy quarks (i.e. 
b= l(UN-2N f ) = 9). 

Let us now write down the spectral representation for the correlator of the trace of the energy-momentum 
tensor: 

<O|T0£(x)C(O)|O> = J 0^p e (k 2 )9(k o )(e- ik *6(x°) + e ik *e(-x )) 



da 2 p e {a 2 )A F {x-a 2 ), (11) 

where the spectral density is defined by 

Pe(k 2 )6(k ) = J2(^) 3 5HPn - fc)IH^|0)| 2 , (12) 

n 

the phase-space integral should be understood in Eq. ([l2|), and 

iA F (x; a 2 ) = i J -^S(k 2 - a 2 )6(k )(e- lkx 6{x°) + e lkx 6(-x°)) (13) 

is the Feynman propagator of a scalar field in the coordinate space. Substituting the representation ( |i"i"| ) 
in Eq. (^), the potential can be expressed as a superposition of Yukawa potentials corresponding to the 
exchange of scalar quanta of mass a: 

V e {R) = -i(d 2 ^) (^-) J dtj da 2 Pe (a 2 )A F (x;a 2 ) 

-(*i)"(x)7^55^- 

Our analysis so far has been completely general; the dynamics enters through the spectral density jl^). 
Let us first evaluate this quantity in perturbation theory, where it is given by the contributions of two-gluon 
states (see Fig. 1(a)) defined by 

pf{q 2 ) = Y,(^) 3 SHpi +P2- q)\(p 1 e 1 a,p 2 e 2 b\6^\0)\ 2 , (15) 

where the phase-space integral is understood, as well as the summations over the polarization (£1,2) and 
color indices (a, b) of the two gluons. The calculation for SU (N) color (see Appendix A) gives 



the appearance of q 4 dependence in Eq. ([16]) is of course natural from dimensional arguments. Performing the 
integration in Eq. ( |l4| ) over the invariant mass a 2 from zero to infinity, we get the following result (N = 3): 

This result can also be derived by the functional method of Bhanot and Pesk in M (see Appendix B). 

Several remarks are in order here: The oc R~ 7 dependence of the potential (17) is a classical result known 
from atomic physics [Q; as is apparent in our derivation, the extra as compared to the Van der Waals 
potential oc R~ 6 is the consequence of the fact that the dipoles fluctuate in time, and the characteristic time 
of fluctuation t ~ e^ 1 (e is quarkonium binding energy) is small compared to the spatial separation of the 
"onia": t<^LR— note an explicit integration over time in Eq. (|l4|). This illustrates, in a somewhat different 
way, the original argument of Voloshin Q] that the physical picture behind the OPE is orthogonal to the 
potential model - the latter is based on the assumption of instantaneous interaction, whereas the former 
is based on the assumption that the internal frequency of heavy quarkonium 1/eo is much higher than 
the frequency of external soft fields. Retardation effects make questionable the possibility to describe the 
interactions of quarks inside a heavy quarkonium by a local potential. In our case, applying the OPE method, 
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(a) (b) 

FIG. 1. Contributions to the potential between quarkonia from (a) two-gluon exchange and (b) correlated two-pion 
exchange. 



we first average the interactions with soft gluons over the quarkonium internal state, which corresponds to 
the infinite retardation. With the resulting coupling between the quarkonium and the gluons, the potential 
description of onium-onium scattering is adequate since at low energies the relative motion of heavy quarkonia 
is slow. The retardation effects manifest themselves in the modification of the shape of the potential. 

We note that although the matrix element of the operator 8£ can in general be non-perturbative, in 
perturbation theory 9£ is of order g 2 , and accordingly the potential ([l7|) has the prefactor g 4 . Then the 

second term g 2 9^ in Eq. (^|), which describes the tensor 2 ++ state of two gluons, gives the contribution in 
the same order in g. Adding this contribution to V$ in Eq. (|l7|), we recover the complete result of Ref. [|| 

" p, < r >-» i ( j 4) 2 5^ ; <i8) 

Note that our d~i is related to the di in Ref. J^] by d2&oeo = <^2<? 2 - This perturbative expression is valid when 
a ,l/e < J R< A Q( t, D . 



III. BEYOND THE PERTURBATION THEORY: THE ROLE OF GOLDSTONE BOSONS 



At large distances, the perturbative description breaks down, because the potential becomes determined 
by the spectral density at small q 2 , where the transverse momenta of the gluons become small. 



A. Broken scale invariance 



To see the importance of non-perturbative effects explicitly, let us consider the correlator of 9% 



U(q 2 ) = i f d 4 xe^(0|T^(x)^(0)|0) = f da 2 } , . (19) 

J J <j z - q l - «e 

An important theorem |nj for this correlator states that as a consequence of the broken scale invariance of 
QCD, 

11(0) = -4<O|0£(O)|O). (20) 

Note that the r.h.s. of Eq. ( po| ) is divergent even in perturbation theory, and should therefore be regularized 
by subtracting the perturbative part. The vacuum expectation value of the 9£ operator then measures the 
energy density of non-perturbative fluctuations in QCD vacuum, and the low-energy theorem (20) implies a 
sum rule for the spectral density: 

J ^-[pf^) - pgV)] = -4<O|0£(O)|O) = -16 £vac ^ 0, (21) 
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where the estimate for the vacuum energy density extracted from the sum rule analysis gives e vac ~ 
— (0.24 GeV) 4 Since the physical spectral density, Pg hys , should approach the perturbative one, pff, at 

high a 2 , the integral in Eq. ( pl| ) can accumulate its value required by the r.h.s. only in the region of relatively 
small a 2 . 

In addition, another sum rule |36] p8| , 

J da 2 pt**(a 2 ) = J da 2 pf{a 2 ) (22) 
is implied by the quark-hadron duality. 



B. Matching onto the chiral theory 

At small invariant masses, the physical spectral density of the correlator (|n]) should be saturated by the 
lightest state allowed in the scalar channel — two pions: 

pV(q 2 ) = E( 27r ) 3(54 (^ +pz - q)KApiMP2)K\o)\ 2 , (23) 

where, just as in Eq. (|l5|), the phase-space integral is understood. 

Since, according to Eq. (poj), 9£ is a gluonic operator, the evaluation of Eq. ( |23] ) requires the knowledge of 
the coupling of gluons to pions. This is a purely non-perturbative problem. Nevertheless it can be rigorously 
solved, as it was shown in Ref. (see also |^|). The idea is the following: at small pion momenta, the 
energy-momentum tensor can be accurately computed using the low-energy chiral Lagrangian, 

£ = L tr d^Ud^ + ~ m 2 jl tr (U + &) , (24) 

where U = exp (2i7r// w ), tt ee ir a T a and T a are the SU(2) generators normalized by tr T a T b = \8 ab . The 
trace of the energy-momentum tensor for this Lagrangian is 



f 2 

6£ = -2 ijj- tr d^Ud^ - m 2 / 2 tr (U + U^) . (25) 

Expanding this expression ( p5|) in powers of the pion field, one obtains, to the lowest order, 

0£ = -d^ a dfTi a + 2mlir a ir a + ■■■, (26) 

and this leads to an elegant result fic[ | in the chiral limit of vanishing pion mass: 

<7r+7r-|^|0)= (7 2 . (27) 

This result for the coupling of the operator 8£ to two pions can be immediately generalized for any (even) 
number of pions using Eq. ( p5[ ) . 

Now that we know the coupling of gluons to the two-pion state, the pion-pair contribution to the spectral 
density (|2^) can be easily computed by performing the simple phase space integration with the result 

^ = si ^ (28) 

in the general case of Nf light flavors, the coefficient 3 in Eq. (p8| ) should be replaced by (N 2 — 1). Again, the 
q 4 dependence comes only from dimensionality. Multi-pion contributions can be evaluated using Eq. (|2q) ; 
we have found that at small invariant masses their influence is small. The dominant contribution at small 
invariant masses a, in which we are primarily interested here, therefore comes from the tttt state. 

Recalling that to the leading order in OPE the scattering amplitude is dominated by the operator ^g 2 ~E a2 , 

we need to evaluate also the matrix element of the second term in Eq. (^|), (0\g 2 6^ \mr) to complete 
our derivation of the scattering amplitude. As we mentioned in the previous section, this tensor operator 
contributes a substantial fraction, 8/23, to the full perturbative result. However, unlike the scalar operator, 

the tensor term is not coupled to the anomaly. The contribution (0\g 2 0^ \irir} therefore is of 0(g 2 ), and 
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FIG. 2. Scalar-isoscalar formfactor of the pion 



is sub-leading in the heavy quark limit. In this limit, we thus come to the following low-energy expression 



(™\±g 2 W*\0) = (^) q -+(H^.n,h. (2-n 



The matrix element in question is therefore known up to a s and m\ corrections. 

The most important correction due to the finite pion mass is the phase space threshold; to take it into 
account, we modify the spectral density in the following way (q 2 > 4m 2 ): 

2 a 2 \ 1/2 

q 1 - Am 2 x 



"■M-aFlVJ 5; (30) 

this expression should be valid at small q 2 . Substituting this spectral density into the general expression 
(0), we get the potential due to the tttt exchange; at large R 

A 2 \ 2 o 1/2 



"""<*> - -( d v) l^J i< 2 -> 4 (iSf* < 31 > 

Note that this potential is not of Yukawa form. The same i?-dependence of tttt exchange at large distances 
was found long time ago by Levy J39| and Klein |l(|. It has been given previously also by Bhanot and 
Peskin but up to an unknown constant. In our approach, the strength of the potential, as well as its 
dependence on the numbers of colors N and flavors Nf - ~ (7V| — l)/(lliV c — 2iV/) 2 , is fixed by the low-energy 
QCD theorems. 

Note also that, unlike the perturbative result ( p"7| ) which is manifestly 0(g 4 ) (besides a factor ((fcoo/eo) 2 ), 
the amplitude (pjh is O(g ) - this "anomalously" strong interaction is the consequence of scale anomalyn. 



C. Dynamical enhancements in the spectral density 

The low-energy theorems |3l]JlC| not only allow us to evaluate explicitly the contribution of uncorrelated 
tttt exchange; they also tell us that this contribution alone is not the complete answer yet. Indeed, the 
numerical analysis shows that the 7T7T spectral density ([[o]) alone cannot saturate the sum rule © - at large 
cr 2 , the physical spectral density approaches the spectral density of perturbation theory, so the integral in 



4 Of course, in the heavy quark limit the amplitude (|3l|) will nevertheless vanish, since ao — » and eo — > oo. 
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FIG. 3. Spectral density of the correlator 
(0|T^(a;)^(0)|0) at low energy (solid line). The mi- 
correlated two-pion contribution is shown in dotted 
line, and the perturbative one with one-loop running 
coupling constant (Aqcd = 200 MeV) in dashed line. 
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FIG. 4. Physical (solid) and perturbative (dashed) 
parts of the integral (^l|) as a function of the upper 
limit, so. The LET states that the difference of the 
two should be equal to the QCD vacuum energy den- 
sity, 16 |e V ac| ^ 0.053 GeV 4 . 



Eq. (|T]) does not get any contribution; at small a 2 , the irir spectral density (|30|), according to the chiral 
and scale symmetries is suppressed by ~ c 4 . The low energy theorems thus require the presence of resonant 
enhancement (s) |3(]] in the ++ rnr, and perhaps multi-pion, KK and r\r\ channels as well. Here we will 
leave the complete multi-channel problem for future investigations, and study only the influence of these 
resonances in the 7T7t channel on the potential between the color dipoles. 



To do this, we define the pion scalar form factor by 
down the spectral density as 



|0fJ|O) = q 2 F(q 2 ) (in the chiral limit) and write 



32tt 2 



4m; 



1/2 



(32) 



It may be illustrative to consider first the idealized case of a sharp a resonance. For simplicity, let us assume 
that the difference between the physical and perturbative spectral densities is due to this a resonance alone, 
and write the spectral density as Pg hys (s) — ^(s) — c 5(s — m 2 ). The LET ( pi] ) then fixes the contribution 
of the narrow a state of mass m a as 



The corresponding potential is of Yukawa type, 

a 2 \2 /4tt 2 



V(R) 



d 2 ^ 



4ttR 



(33) 



(34) 



In this idealized situation, the strength of the potential is directly related to the energy density of non- 
perturbative QCD vacuum. Note, however, that this simplified model of the sharp er resonance is inconsistent 
with the asymptotics derived from the broken chiral symmetry (Cf. Eq. (|3l|)). 

The formfactor F(s) is directly related to the experimental irir phase shifts by the Omnes-Muskhelishvili 
equation 41 42]. Within the single-channel treatment F(s) has a solution, 



F(s) — exp 



ds' . , 

4m 2 s'{s' 



(35) 



where Sq(s) is the phase shift of the 7T7t scattering in the scalar-isoscalar channel, and formally s\ 
With this formula we can make a full use of the experimental information on the 7T7r correlations. 
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FIG. 5. Integral of the duality relation. The notations are the same as in Fig. ^| 



In our calculation we use a simple analytic form p3[ for the phase shift Sq(s) which has been shown 
to fit the experimental data up to s^tt — 1 GeV 2 . Beyond this energy one should take into account the 
contributions of other channels, such as KK . We performed the integral in Eq. ( |35| ) numerically up to 
Si=(5 GeV) 2 by extrapolating the low-energy fit of the phase shift. When we change s\ to (20 GeV) 2 , the 
change in F(s) at 1 GeV 2 is a few percent. In Fig. ^ we show the resulting scalar formfactor of the pion, 
F(s). The structure of F(s) may be interpreted as due to a broad a and narrow /o resonances. For a 
more realistic evaluation of the formfactor, the multi-channel calculation has to be done; the results will be 
reported elsewhere. 

In this paper, as a simple model for the Pg hys , we will take the form 



Pe ( )_i P?(S) (80 < 6), 



where Sq is a matching scale. 



D. The analysis of the sum rule 

Let us consider the sum rule ( plf ) within our simple model for the spectral density. When the model 
( |36| ) for Pg hys {s) is used, the upper limit of the integral in Eq. ( |2l| ) can be replaced by sq. In Fig. || we 
show the physical and perturbative parts of the integrand in the sum rule ( pl| ) with solid and dashed lines, 
respectively. Since for the spectral density in the perturbation theory there is no scale other than s, the 
coupling constant should be taken running with this scale: 



where a s {s) = An/(b ln(s/A^ CD )) with A QCD =200 MeV. 

We note that the spectral density for uncorrelated pions pg 77 (|3(]) , which is shown in dotted line in Fig. ||, 
has the same functional form as the pg , up to the logarithm and the threshold factor. As a consequence 
of this, we find that the uncorrelated 7T7T spectral density (pOj) cannot obey both the sum rule p|) and the 



duality constraint (22) 



On the other hand, the spectral density obtained with the Omncs Muskhelishvili solution has a non-trivial 
structure (see Fig. ^) ; one can clearly see a narrow peak at the /o resonance region with a shoulder coming 
from the broad "a" around 0.6 GeV. 

We plot the integral in the LET for the physical (solid) and perturbative (dashed) parts separately as 
a function of the upper limit, sq. One can see that the value of the integral for the physical spectral 
density is mainly accumulated in 1 GeV region, and the "cr" contributes to it about 20 % (see Fig. |). The 
perturbative part behaves as s 2 up to the logarithm, weighting the higher energy region. The LET tells us 
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FIG. 6. Potential ( |l4| ) between two J/ip's (solid line). Contributions of the spectral densities of s > so and s < so 
with y^so = 2 GeV, respectively, are plotted in dotted and dashed lines. The perturbative result (^) is shown in 
dashed-dotted line, for reference. 



that the difference of these two contributions should be equal to the energy density of the QCD vacuum. 
In our model for the spectral density, the LET (|2l]) is satisfied when we choose sq=(2~2.5 GeV) 2 . As for 
the duality relation (|2^) , the equality of the integrals of the physical and perturbative spectral densities is 
achieved when we choose so ~ (2 GeV) 2 (Fig. fig5) - this value of the matching scale therefore provides a 
consistent solution to both the LET and the duality relation. 

Even though our spectral density cannot be taken seriously in the high mass region beyond ~1 GeV, our 
calculation nevertheless shows the following: Non-perturbative dynamics of QCD generates enhancements in 
the intermediate mass re gion in the form of hadronic resonances, which make the physical spectral density 
consistent with the LET (21). The narrow /o(980) is more important for the LET than the low mass, broad 
a resonance. Therefore, to discuss the influence of heavier resonances (like /o(1500)) we need to perform a 
coupled-channel analysis including the KK and other states. In the rest of this paper we put the matching 
scale yfso = 2 GeV. 



IV. THE POTENTIAL BETWEEN COLOR DIPOLES 



As a concrete example, let us consider the potential between two J/ip's at rest. Although the charm quark 
is perhaps not heavy enough to justify the heavy quark limit, we try to extrapolate our result to J /if) and 
discuss its implications. 

For the pure Coulombic bound state, e?2 = 7/36, and the Bohr radius and Rydberg energy are given by 
ciq = 4/(3a s m) and eo = (3a s /4) 2 m = l/(agm), respectively (a s = g 2 /4Tr). We have a s (J/^)=0.87 and 
ao=0.20 fm for the J/tp with the phenomenological inputs, eo = 2A/d — M( J/-0)=642 MeV and m=1.5 GeV. 
These values show the application to J/tp will be qualitative at most, because of the large a s value and 

because ao > s ; the latter means that nonpcrturbative effects penetrate inside the radius of J/tp. 

In Fig. H the resulting potential ( |l4] ) between two J/ip's is shown as a solid line. In our model for the 
spectral density, the potential consists of two components, high q 2 (dotted) and low q 2 (dashed), separated 
by So, which we set (2 GeV) 2 . (As in the heavy quark limit, we omit the contribution from the tensor 

exchange, # cP)- 

First, we see that the potential at large distances is naturally determined by the spectral density of the low 
q 2 region. Moreover the total strength of the potential at large distances is enhanced by the non-perturbative 
spectrum of QCD, compared to the formal perturbative result ( |l7|) denoted by the dashed-dotted line. The 
region where the two components compete is R ~ 0.5 ~ 0.6 fm, which is much larger than the scale 

I/O , . 1/2 

determined by s ~ 0.1 fm. This is in contrast with naive expectation that beyond the scale s , the 
potential should be dominated by the non-perturbative spectral density. The reason for this lies in the large 
value of a s (J/ip), reflecting the fact that charmonium is still far from the heavy quark limit. 

In the discussion of the LET (pin, we used the running coupling constant, while the coupling constant 
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used here is frozen at the J ftp scale. This is because in the heavy quark limit, it is natural to renormalize the 

coupling constant at the scale of quarkonium, <?(eo), with eo 3> \ Aqcd- The matrix element of G 2 should 
then contain the effects of quantum fluctuations below this energy scale. Within the perturbative approach 
the renormalization group ensures independence of the final result on the choice of renormalization point, at 
least in the leading- log approximation, which we used in Eq. (^). In the case of J/ip, the renormalization 
scale (chosen at the binding energy, eo) is still lower than So, and the spectral density ( |l6| ) with fixed a a {J/ip) 
is significantly larger than the one with the running coupling constant (|37]). Again, this reflects the fact 
that non-perturbative effects penetrate "inside" the J/tp- The most important feature seen in Fig. |g| is the 
dominance of low-q 2 enhancements in the spectral density in the behavior of potential at large distances. 



V. THE SUM RULE FOR THE POTENTIAL 



We can derive an interesting sum rule for the strength of the potential; according to Eqs. (|14|) and (|17|), 
we have 

jy^{V e (R)~Vr(R))=-{d 2 a l)\^) 2 J d -^{ p fy^)- P f{a 2 ))T{2,aa) , (38) 

where a should be chosen to be of the order of the onium radius, and T(z 7 p) — dtt z ~ 1 e~ t . As we discussed 

previously, the physical spectral density p^ ys {& 2 ) differs from the perturbative one, pf'io 2 ), in the region 
a 1 < so- In the heavy quark limit, a cx l/(a s m) and a^/io <C 1. Therefore we can re-write the sum rule 
( |38| ) in a more suggestive form: 

£°*b.(v 9 w-v?(r)) =-(^f) 2 ( ! r)7 ^(/>f y V)-P?V)) 

= -(d" 2 £) 2 (^) 2 16| evac |, (39) 

which relates the overall strength of the interaction between small color dipoles to the energy density of the 
non-perturbative QCD vacuum. 



VI. QUARKONIUM INTERACTIONS WITH PIONS 



As another application of our formalism, we evaluate the cross sections of elastic scattering 7r<E> — > 7r$ and 
of excitation process 7r$ — > 7r$'; the latter cross section was previously computed in Refs. pO,^T|. These 
cross sections are important for the analyses of quarkonium production in heavy ion collisions [32^33|. The 
fact that soft pions effectively decouple from heavy quarkonia was previously noted in Rcf |19[| . 



A. Elastic 7r$ scattering 



Within the OPE formalism it is straightforward to write down the amplitude of pion-quarkonium 
elastic scattering at small energies: to the leading order in OPE, 

M kl (P', P ';P,p) = -d 2 ^(n k (p')\lg 2 E a2 (0)\n l (p)). (40) 

The matrix element (7r fc | 5 2 E a2 (0)|7r z ) can be found from (ttV | 5 2 E a2 (0)|0), (§|), by crossing; the LET (||) 
tells us that up to a s and m 2 , corrections 

(^V)l^ 2 E a2 (0)|7r'b)) = Wk'(p)) = 6 kl ^f t F(t), (41) 

where t = [p — p') 2 , and we have introduced, as before, the pion scalar formfactor, F(t). Taking into account 
the non-relativistic normalization of the $ state, we have the expression for the total elastic cross section in 
the CM frame, 
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FIG. 7. nj/ip elastic cross section (solid) as a function FIG. 8. TrJ/ip —> tti/j' cross section. The notations are 

of the CM momentum. Dashed line is the case of F(t)=l. the same as in Fig. Fl 



ZpoVrel J (Z^r J ( 2?r ) 2 Po 

where v re i — y/ (P ■ p) 2 — M 2 m 2 / P°p° is the relative velocity of the incoming J/tp and pion, and p 2 = 
(s - (M - m n f)(s - (M + m 7r ) 2 )/4s is the CM momentum. 

The result for the elastic nJ/'ip cross section is shown in Fig. [t]. The pion scalar formfactor F(q 2 ) and 
other parameters are the same as in the previous section. Note that the chiral symmetry requires a strong 
momentum dependence, t 2 in Eq. (|42|); therefore at low energies the J/ip interaction with pions is very weak. 
Extrapolation of the scalar formfactor F(t) to the scattering region, t < 0, induces additional suppression 
of the tt J/tp interaction. At small energies (see Fig. @) the cross section is on the order of 0.01 mb, which 
is much smaller than the geometrical cross section of the J/i/j. For quarkonium production in heavy ion 
collisions, this implies that the interactions with secondary pions do not contribute to the broadening of the 
quarkonium transverse momentum spectra. 



B. 7r«I> — > 7r$' transition amplitude 

Our next example is the transition process, 7r$ — ► 7r<£>'. In this case, however, the transferred momentum 
is on the order of the binding energy, A = M' — M = O(eo), which may invalidate our assumption on 
the factorization between the short and long distances. Fortunately, since the size of quarkonium ao ~ 
l/(g 2 m) <C 1/eo ~ l/(g 4 m) in the heavy quark limit, the double-dipole form 



M = (0V|tr 



r • E -r ■ E 

H a +e + iD° 



l<M (43) 



is still valid S. The structure of this formula is transparent: the initial quarkonium <!> absorbs/emits a gluon, 
then propagates with the internal energy, — e + Q, and emits/absorbs another gluon to form a color-singlet, 
excited quarkonium state these gluons originate from pions. 

To apply our formalism, let us approximate —iD° in Eq. ( |43| ) by the typical value of the gluon momentum, 
A. Within this (rather crude) approximation, the quarkonium part and the pion part can be factorized in 
the matrix element, (p3|), and the relevant Wilson coefficient, which for this process reads 

4- - ^rM'V 1 - r^). (44) 
£o oN H a + e ~ A 

The transition matrix element then reduces to the same form as in the elastic case: 
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FIG. 9. Schematic picture of the potential between heavy onia. 



4n* 



(45) 



and the total transition cross section can be written as 



<T(S 



1 MM' 

167TS p 2 



3* 



4^2 



dirt) t 2 \F(t)\< 



(46) 



In Fig. || we show the result for the n J/ ip — ► nip' cross section, evaluated assuming the Is and 2s Coulomb 
wave functions, and A = (3/4)eo- It shows that the cross section is on the order of 0.01-0.1 mb. 

We also evaluated the partial width of the ip' due to ip' — > ipnn decay within the same formalism, and 
obtained T=260 (70) keV with (without) using the formfactor F(s). This should be compared with the 
experimental value, 135 ± 20 keV J44| . We conclude that our calculations, due to the assumption of the 
heavy quark limit, hold to within a factor of 2 only. Additional confirmation of the s dependence of the 
matrix element comes from the dipion invariant mass distribution in the ip' — > ipnn decay J45[ . 



VII. SUMMARY AND DISCUSSION 



We have shown that at large distances the interaction of small color dipoles becomes totally non- 
perturbative. This result has a deep physical origin: indeed, one can trace it back to the sum rule @ 
for the correlator of the energy-momentum tensor, which reflects the fact that the non-perturbative vacuum 
of QCD is characterized by non-zero energy density. 

For QCD practitioners, "non-perturbative" is often a substitute for "incalculable". Nonetheless, in our 
case, we were able to evaluate explicitly this, non-perturbative, scattering amplitude in a model-independent 
way. The key ingredients in our approach were i) the use of spectral representation in the t— channel and 
it) the low-energy theorem arising from the (broken) scale and chiral invariances of QCD. 

What are the implications of our results? First, we find that the long-distance interactions of small color 
dipoles are dominated by pion clouds; the qualitative picture of this interaction is illustrated in Fig. ^. 
The size occupied by the heavy quark-antiquark pair in the quarkonium (see Fig. ^ is a ~ l/(g 2 m); 
the gluon cloud spreads up to the distances on the order of the inverse binding energy 1/eo ~ l/(g 4 rn), 
since the typical momentum K of gluons is K ~ eo. (This picture of quarkonium structure emerges also 
from the NRQCD approach of Bodwin, Braaten and Lepage p^] .) The pion cloud begins to dominate at 

the distances ~ s Q , and spreads up to the distances (2/i) -1 , where fi is the pion mass (sq is the mass 
scale at which the non-perturbative effects begin to dominate in the spectral density, see Section III). This 



pion cloud may as well be important in high-energy scattering. One may even speculate (see Bjorken 28 



that pions are responsible for the so-called "soft pomeron" , even though it is not yet clear how to extend 
our calculations to high-energy scattering - this would require evaluation of higher orders in the multipole 
expansion. Nevertheless the possibility that the diffusion of partons toward small k t makes pionic degrees of 
freedom important looks very plausible to us. 
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Second, the fact that pions (and therefore light quarks) dominate the long-distance interactions of heavy 
quark systems is important for lattice QCD simulations. Even though naively one may think that light 
quarks are relatively unimportant for the studies of heavy quarkonia on the lattice, our findings show that 
the opposite is true. This suggests that to extract the properties of heavy quarkonia from the lattice QCD 
one has to use "unquenched" theory with light quarks. The importance of pionic degrees of freedom in 
determination of the mass splittings of heavy quarkonia was investigated in Ref . 

Third, we find that both inelastic and elastic TrJ/ip scattering cross sections are very small, less than 0.1 
mb. The smallness of inelastic cross section suggests that pions are very ineffective in dissociating J/ip's, 
lending support to the idea to use quarkonia as a signal of deconfinement |3|] . The smallness of the elastic 
cross section explains why the transverse momentum distributions of J/t/j's seem to be unaffected by the 
final state interactions with secondary pions, whereas much bigger tp n s, to which our multipole expansion 
analysis does not apply, can show significantly larger mean transverse momenta. 
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APPENDIX A: DERIVATION OF EQS. (14) AND (16) 

The Feynman propagator of a scalar field cp of mass a is defined by 
iA F {x;a 2 ) = i(T<p(x)<p(0)) 

'' k S(k 2 - <T 2 )8(k )(e- tkx 8(x Q ) + e lkx 8(-x°)) 



(27T) 

d 4 k e- lkx 
(2tt) 4 a 2 - k 2 - ie 



f K 1 (^/(-x 2 +ie)a 2 ) , (Al) 



4tt 2 y/(- x 2 + i e )a 2 

where K\ is the modified Bessel function. The Born amplitude of one-boson exchange with coupling g is 

iM B (q) = / d 4 xe tcl - x (Tig<p(x)igtp(0)) = ig 2 [ d 4 xe^iA F (x; a 2 ), (A2) 



which may be related to an interaction potential by iA4(q) = —iV(q). Going to Euclidean space, where 
x\ — t 2 + x 2 = t 2 + i? 2 , we can show that 

/oo />oo -i 2 i 2 

dtiA F (x;a 2 ) = ~g 2 / dr — ~^=K, (Jx 2 E a 2 ) = -/^e"' fl . (A3) 



To the leading in the OPE, the potential between color dipoles (14) is a superposition of this Yukawa 
potential with the spectral function, pg(a 2 ). In the perturbative calculation, the matrix element of G 2 
between the vacuum and the two-gluon state is 

( Pl e 1 a, P2 e 2 b\G a0c G c ap \O) = 4(-p x ■ p 2 e\ ■ e* 2 + Pl ■ e* 2 p 2 ■ el)5 ah + 0(g 2 ), (A4) 

where pi, e, and a are the momentum, polarization and SU(N) color index of the gluon, respectively. Noting 
that the sum over physical polarizations yields the projection, J2 po i £ n £ t — ~PiPj/P 2 : we have a compact 
expression, 

\{ P ieia, P2 e 2 b\G 2 \0)\ 2 =4 2 (N 2 -1) 2( Pl ■ p 2 ) 2 = 8(N 2 - l)q 4 , (A5) 

pol,col 

which depends only on q 2 — (pi + p 2 ) 2 . With the phase space factor of two identical particles, 
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1 f d 3 Pl f <l A p, 



2 J (2.)3 2wi y (2 7 r)32. 2 ^> 1+ ^-^ = 32^' (A6) 

we find Eq. ( |l6|) as the spectral density of the correlator of 8£ = — (bg 2 /32TT 2 )G 2 , to the leading in g. 

Similarly the spectral density of the two-pion states ( p8| ) can be calculated, but with (7T fc 7r'|#^|0) = q 2 S kl 
in the chiral limit (k, I = 1, 2, 3). 



APPENDIX B: ALTERNATIVE DERIVATION OF EQ. (17) 

To confirm our result, we derive here Eq. ([l7]) applying the functional method of Ref. || to the scalar 
part of the interaction (Q) (in Euclidean space); 

7.? v 2 r°° 



Vo(R)=-(d 2 ^) J J T (0\^G 2 (x)^G 2 (0)\0) 



/_ 2 \2„4 poo 

= "(*^) §5 / ^t(0|G^(*)G^,(0)|0} 2 , (Bl) 

where all indices are summed over. Using the expression for the two-point function of the gluon in Feynman 
gauge, 

(W A l(y)) = ^j^y^^ ab > (B2) 



we have 



(G^i/ i/' (0)) = 4^2"^ {<W(<W a;2 ~ ^ x n x w) + <W(<W£ 2 - 4x v x V ') 

-Sv^iS^x 2 - Ax^Xv,) - 6^(5 Ufl ,x 2 - 4x u x^)} +0(g 2 ), (B3) 

and 

24 1 

7T 4 X 8 



(G^)G^,(0)) 2 = (iV 2 0( 5 2 ). (B4) 



Substituting Eq. ( |B4[ ) into Eq. (Bl), we obtain the leading expression ( |l7| ) for the potential of the scalar 
part (x 2 = t 2 + R% 
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